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Abstract: A separator failure during a store in 2002 led to a drop in luminosity, to in- 
creased emittance growth and to a drop in beam lifetimes. We show that a simple diffusion 
model can be used to explain the changes in emittance growth and beam lifetimes. 
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1. Introduction 

Emittance growth of beams when they are in collision occurs due to many sources: beam- 
beam interactions, magnetic nonlinearities, intra-beam scattering, scattering off the resid- 
ual gas and possibly others. The dynamics of the emittance growth is complicated and it 
depends strongly on the tunes. It is not always clear that the dynamics can be described by 
a diffusion process at all particle amplitudes in each beam. However in one store early in 
Run II, there was a sudden drop in a separator voltage in the Tevatron and the subsequent 
enhanced emittance growth and intensity lifetime drop could be described by a simple 
diffusion model. In this report we analyze the luminosity drop, compare the measured 
value with the expected drop and analyze the change in beam lifetimes. We show how a 
simple model of diffusive emittance growth and a change in physical aperture provides a 
quantitative explanation for the change in lifetimes. 

2. Separator failure and luminosity drop 

After about 13.5 hours into store 1253 (April 26, 2002), the voltage on the bottom plate of 
the horizontal separator at A49 dropped from a value of -90kV to -25kV. This immediately 
lowered the luminosity at CDF and DO. The emittance growth rate increased and the 
lifetimes of both protons and anti-protons fell. Table [I] shows some of the key beam and 
machine parameters. 

Figure [1] shows the measured luminosity drop at the CDF experiment after the separa- 
tor failure. The luminosity at DO dropped by a similar amount. 
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P* at CDF, DO 




0.35m 


Horizontal tune 




Vx = 20.585 


Observed Luminosity drop at CDF 




41.4% 


Observed Luminosity drop at DO 




42.3% 


Total proton beam intensity before drop 




5.78 xlO^^ 


Initial emittances (p and pbar) £x, £y [;rmm-mrad] 




22,21 


Estimated final emittances {p,p) [TTmm-mrad] 




(26 - 30, 25 - 29) 


Average emitt. growth rate [;rmm-mrad/hr] 




0.3-0.5 


Length of store [hrs] 




15 


Location of BPMs nearest to CDF [m] 


7.5 


upstream and downstream 


BPM upstream of CDF 


^= 


159.5 m, \j/x = 27tx 20.337 


BPM downstream of CDF 


/3x= 


160.44 m, Yx = 2nx 0.238 


A49H Separator 




867.67m, Xj/^ = 2nx 20.329 



Table 1. Relevant beam and machine parameters in Store 1253. 
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Figure 1. Average CDF Luminosity (raw data) and fits before and after the separator failure. The 
luminosity at CDF dropped by 41% right after the separator failure. 



2.1 Closed orbit shift at CDF and DO 

Here we estimate the orbit change as a consequence of the change in separator voltage at 
A49. The shift in the closed orbit due to a kick AO can be found from 

^ , ^ _^2^f£P!^cos[nv,-Ms) - Y^epl] (2.1) 
2 sm KVx 

where fisep, j8 (s) are the beta functions at the separator and at s respectively, \if{s) , Xifsep are 
the phase advances from a reference point to s and the separator respectively and is the 
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horizontal tune. The kick resulting from the electric field S across the separator plates of 
length L on a particle of energy E is given approximately by A0 = (SL/E. At 980 GeV, a 
change in voltage of 65 kV across the separator plates with a gap of 5cm results in a kick 
of 

A0 =3.41/irad (2.2) 

Using the above expressions we find that the proton's horizontal closed orbit would 
move by 

AxcoiCDF) \sep = 29.9Mm, Axco{DO) \sep = -30.3Mm (2.3) 

The orbit separation dc^ between the beams at the IPs would be twice the above value if 
the protons and anti-protons undergo the same but opposite changes in orbit. However the 
beam-beam kick with separated beams also induces an orbit shift and this will be larger 
for the anti-protons. Calculation of the beam-beam induced orbit kick requires that we 
know the separation between the beams, but that is precisely the quantity that we want 
to predict. We will approximate the beam-beam induced orbit kick by assuming that the 
beam separation is twice the shift in the proton orbit. In that case, the beam-beam kick 
assuming round Gaussian beams is 

= ^^^^ (. + j)2+,2 { 1 - -Pt-^-^^] I (2.4) 

where ^ is the beam-beam parameter, e is the un-normalized rms emittance and d is the 
horizontal separation between the beams. Extracting the dipole part of the kick, 

HmbmM = 8;r^el |l -exp[-^]| (2.5) 

Since the sign of the pbar orbit offsets at CDF and DO due to the separators have opposite 
signs, the kicks experienced by the anti-protons due to the dipole beam-beam kicks at CDF 
and DO have opposite signs. Hence the contribution of the beam-beam kicks at CDF and 
DO to the orbit shift at CDF is 

AxcoiCDF) \bb = ^ '. [cos nv, - cos{nv, - \ v/,(DO) - V^x(CDF) |)] (2.6) 

2 sm nVx 

We have assumed here that the beta functions at the IPs did not change much. 

From the average proton bunch intensity ofA'p = 1.61xlO^^ and an expected proton 
emittance of 30;r at this stage of the store (this number is found later by a self-consistent 
calculation), we find that the beam-beam parameter for anti-protons at this stage was 
^ = 3.92 X 10 Hence the beam-beam kick using the value of the orbit offset found 



in Equation ( |2.3[ ) is 

AxU,„(CDF) = 5.14Mrad (2.7) 

while at DO, the kick has the opposite sign. Note that this kick is larger than the kick due 
to the change in the separator voltage, cf. Eq. (|2.2|). However because of the small beta 
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function at the IPs, the change in orbit due to these beam-beam kicks is quite small, 



Axco{CDF)\bb = -0A2^m (2.8) 

using Equation (|2.6|). This is almost two orders of magnitude smaller than the orbit shift 
due to the separator and can be neglected. 

The predicted luminosity in terms of the luminosity .5fo before the separator failure is 
found from 

where dco ~ '2Axco\sep- This calculation depends on the emittances at the time of the 
failure. With the initial emittances and average emittance growth rates shown in Table 
1, the proton emittances after 13.5 hrs were likely to be in the range (26, 30);rmm-mrad 
and anti-proton emittances (25, 29);rmm-mrad. We find using the orbit shifts in Eq. ( |2.3D 
(virtually the same at CDF and DO) and the low and high end of the emittance range, 

= 0.47, {£p = 26, Epbar = 25) 

^ = 0.42 {£p = 30, £pbar = 29) (2.10) 

These values are to be compared with the observed relative drops in luminosity of 0.414 
at CDF and 0.423 at DO. This suggests that the emittances were more likely at the higher 
end of the quoted range. 

Another test of the optics is to propagate the measured orbit changes at the BPMs 
closest to CDF back to CDF using 



j5{s2) cos[;rv - 


IVisi) - Vsep\] 


j5{si) cos[;rv - 


-¥sep\] 



Ax{s2) = " i M^i) (2.11) 



where si is the location of a BPM and S2 is the location of an IP. This expression does not 
depend on the kick angle at the separator nor upon the beta function at the separator. 

From Figure ^ we observe that just before the failure, the proton horizontal position 
was relatively steady at 1.349 mm, then it falls for about 15 minutes after which it stabi- 
lized at 1.217 mm. Similarly at the downstream BPM, the horizontal positions at these 
same times are -6.375 mm and -6.214 mm. This slow decay in the position is related to 
the long integration time (about 15 mins) of these collision point monitors(CPMs). 

The observed shifts at the CPMs around CDF were 

Upstream : Ax'y' = -0. 132 mm Downstream : Ax^^' = 0.161 mm (2.12) 

Propagating these orbit shifts to the IP, we find using the upstream BPM that the expected 
orbit shift at CDF and relative luminosity drops are 

Aj^ 

AxcoiCDF) = 0.0285 mm, -^{CDF) = 0.384 (2.13) 
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Figure 2. Proton horizontal beam position at BPMs upstream and downstream of CDF. The jump 
in position coincides with the separator failure 



while using the downstream BPM, we find 
AxcoiCDF) = 0.0241 mm, 



{CDF) = 0.293 



(2.14) 



This calculation shows that the upstream BPM was more consistent with the observed 
luminosity drop. This may indicate either more errors downstream from the IP to the 
CPM (this is unlikely since there is only the detector between the IP to the CPM) or that 
this downstream CPM reading was less reliable. 



3. Lifetimes and emittance growth times 

The luminosity in terms of beam parameters is 

where frev is the revolution frequency, is the number of bunches in each beam, Np.Np 
are the proton and anti -proton bunch intensities respectively, eN,p^ £n,p are the 95% emit- 
tances of the beams, Os is the rms bunch length and is the hourglass form factor 

Jf{z) = Vnze'\l-^{z)), 4>(z) = ^ fe-'^dt (3.2) 

JO 

The average bunch length in Store 1253 was 2.6 nano-seconds or a, = 78cm. With j8* = 
35cm, the hourglass reduction factor for z = P* /Os = 0.45 is J^{z) = 0.51. 

The luminosity lifetime can be calculated from the beam parameters by taking the 
logarithmic derivatives. Defining the luminosity and intensity lifetimes as 

1 _ I 1 _ 1 dNp 1 _ 1 dNp 

^ dt ' Tp Np dt ' Xp Np dt 
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Figure 3. Proton bunch intensity (left) and anti-proton bunch intensity (right) before and after the 
separator failure. Note the change in slope at around 13.5 hours 



while the longitudinal bunch length and transverse emittance growth times are 
1 1 dOs 1 Id 



(Tj dt 

Then the luminosity lifetime is 

1 



£n,p + £n,p dt 



i£N,p + £n,p) 



(3.4) 



z dJ^ 1111 
+—+—+— 



M' dz Is tp 



(3.5) 



Here z = /Os- Using the expression for the hour-glass form factor J^, this can be 
rewritten as an expression for the emittance growth rate, 



1 



1 



(1 



1 1 



(3.6) 



We calculate the emittance growth time Tg^ from the measured values of the other time 
scales. This is useful because at this stage in Run II the synchrotron light monitor was 
not available, so there was no direct measurement of the transverse emittance growth rate 
during the store. 

Figure ^ shows the proton and anti-proton bunch intensities an hour before and after 
the separator failure. There is a clear change in the intensity lifetimes before and after the 
failure. 

Figure ^ shows the bunch length as a function of time over the store. There is evident 
growth in the bunch length over the 15 hours of the store but from the data then available 
(bunch lengths every 15 minutes) it was not possible to discern a change in the growth 
of the bunch length after the separator failure. We assume that the growth rates of the 
bunch length were the same an hour before and an hour after the failure. Tabled shows the 
measured lifetimes and growth times and the calculated transverse emittance growth time 
before and after the failure. 
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Figure 4. Proton bunch length during the store. The separator failure occurred about at 13.5 hours 
into the store. 





Before separator failure 


After separator failure 


Luminosity lifetime [hrs] 


17 


10 


Proton lifetime [hrs] 


198 


100 


Anti-proton lifetime [hrs] 


64 


36 


Average Bunch length [nsec] 


2.6 


2.6 


Bunch length growth time [hrs] 


81.5 


81.5 


Transverse Emittance growth time [hrs] 


24.5 


15.3 



Table 2. Lifetimes and growth times from data 1 hour before and an hour after the separator failure. 
The transverse emittance growth rates are calculated from Equation (3.6). 

3.1 Beam lifetime 

The separator failure changed the beam orbits around the ring. The changes in the proton 
orbit were calculated using a MAD optics file. The maximum horizontal orbit change was 
1.8(7 radially outwards while the rms orbit change was 0.59a. In the vertical plane, the 
maximum orbit shift was about 1.8(7 downwards but the rms orbit change was smaller, 
0.24a, as expected. An important consequence of the change in orbit is a change in phys- 
ical aperture. Since the collimators define the limiting aperture, movement towards them 
would reduce the physical aperture. Figure |] shows the simulated proton horizontal orbit 
change at several collimators in the ring. The largest change at a collimator location was 
about 0.73a towards the F17 collimator. 

We have seen that the emittance growth rate increased and the orbits changed signif- 
icantly after the separator failure. Could these two phenomena explain the sharp drop in 
beam lifetime? 

We assume that the beam density distribution function evolves according to the dif- 
fusion equation. For simplicity we will consider the transverse distribution function can 
be decoupled as the product of horizontal px and vertical py distribution functions and 
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Figure 5. Change in the horizontal proton orbit (in units of the rms beam size) at the collimators 
due to the separator failure. The largest change was at the F17 collimators where the protons moved 
by about 0.73 a to the radial outside. 



consider only the evolution of Px, 



dt 



d 

dTx 



dpx 



(3.7) 



Here D{Jx) is the diffusion coefficient in the action = [x^ + (Pxx' + ax)^]/l5x. The lim- 
iting physical aperture is assumed to be in the horizontal plane at an amplitude and 
corresponding action at the aperture is Ja =A^/Px- Under the assumption of indepen- 
dence of the transverse planes, the total number of particles in the beam at time t can be 
written as N{t) = NoNx{t)Ny{t) where A'o is the initial number of particles, A'^:^ is a scaled 
time dependent number defined by 



Ja 



Px{Jxj)dJx, Nx{0) = l 



(3.8) 



and a similar expression for Ny. We assume that particles that diffuse out to the aperture 
Ja are lost. 

The beam emittance £x is related to the average action which is defined as 

{Jx)= f^^ JxPciJxj) dJx (3.9) 
Jo 

where pc is the conditional density which accounts for the particle number changing in 
time and hence is defined as pc = Px/Nx- From the definition it is clear that pc is nor- 
malized to unity at all times. It follows from the diffusion equation that the average action 
evolves as 

-(7,)=/ D'{Jx)pcdJx+[JA~{J.)]D{JA)^{JA) (3.10) 
at Jo oJx 
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If the density falls sufficiently rapidly and the aperture is far enough away from the beam, 
then dpc{JA) / dJx — ?■ and the second term in the above equation can be dropped. With 
this simplification, 

^{J.)= ['D'{J,)pcdJ, (3.11) 
at Jo 

If the diffusion coefficient increases linearly with the action, 

D{J,)=DoJ, (3.12) 

then it follows 

tM-)=Do (3.13) 
at 

The effective transverse emittance calculated from the luminosity was found to increase 
nearly linearly with time during stores, we may therefore assume that Equation ( |3.12| ) is 
valid, at least in the core of the beam. 

We can now calculate time scales associated with the diffusive motion. The mean 
escape time for particles to travel to the absorbing boundary at Ja is defined as [Q] 

dJ, (3.14) 



Assuming Equation ( |3.12D , we obtain 



D{J, 



tesc=^ (3.15) 



This escape time is related to the beam lifetime. 

The lifetime can be calculated by a more complete analysis as in Edwards and Syphers 
[jlp. The diffusion equation can be solved analytically when the diffusion coefficient is 
linear or quadratic in the action. With a linear dependence as assumed above, the density 
at time t is 



pAJx,t) = Y.CnMX,J-^)exp[-^^] (3.16) 

n \ JA ^ JA 

where 

^« = Tin^j rVo(^)^o(A„./f ) dJ, (3.17) 

Jf{An)JA JO J A V J A 

Jo, Ji are the zeroth and first order Bessel functions, the A,/s are the n'th roots of 7o and po 
is the initial density. For an initially Gaussian distribution in phase space, the distribution 
in action is an exponential, 

ccJ 

Po(7x) = aexp[ — , a = (3.18) 
Ja 2a^ 

Assuming that the beam is sufficiently far from the aperture, the coefficients c„ simplify in 
this case to 
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Keeping only the first and dominant term in the solution for the density Equation ( |3.16| ), 
the scaled partial number of particles A^^ in the beam simplifies to 

We assume that the lifetime was determined by particles reaching the horizontal aperture 
first. It follows that the lifetime defined as 

tL = — = — ~0.7 — (3.21) 

dN^/dt XI Dq Do 



This expression is very close to the mean escape time tesc calculated in Equation ( [3.15D . 

We will now use Equation ( ]3.21| ) to relate the beam lifetimes before and after the 
separator failure. Equating 

Z)o - ^ = ^ (3.22) 

dt Te 

where Cq is the initial emittance and Tg is the emittance growth time. Then 

fL = 0.7:^Te (3.23) 

Before the separator failure, the beam aperture was approximately 6a at one or more of the 
collimators. Hence before the failure, Ja = (6a)^//3^. After the failure the beam moved 
closer to the physical aperture by 0.7a. Hence Ja = (5.3a)^/ after the failure. 

^ (after) ^ 7a (after) Te (after) 
?L(before) 7a (before) Te (before) 

From Table we find that Te (before) = 24.5hrs and Te (after) = 15.3hrs. Hence 

A^^ = (^)2l^=0 49 (3 25) 

rL(before) ^ 6 ^ 24.5 ^ ^ ^ 

From Table we find that the ratio of the measured lifetimes is = 100/198 = 0.51. Hence 
the predictions of the one dimensional theory are in very good agreement with the observed 
drop in lifetime. The increased emittance growth after the separator failure may have been 
due to multiple sources resulting from the change in orbits including change in tunes, 
increased long-range beam-beam effects from smaller separations at some locations and 
larger nonlinear fields in some magnets etc. 

No matter what the sources of emittance growth were, we have shown that the drop 
in lifetime after the separator failure was consistent with a simple model of diffusive emit- 
tance growth and the beam center moving closer to a physical aperture. 
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